First-principles phase-coherent transport in metaUic nanotubes with reahstic contacts 
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We present first-principles calculations of phase coherent electron transport in a carbon nanotube 
(CNT) with realistic contacts. We focus on the zero-bias response of open metallic CNT's considering 
two archetypal contact geometries (end and side) and three commonly used metals as electrodes (Al, 
Au, and Ti). Our ab-initio electrical transport calculations make, for the first time, quantitative 
predictions on the contact transparency and the transport properties of finite metallic CNT's. Al and 
Au turn out to make poor contacts while Ti is the best option of the three. Additional information 
on the CNT band mixing at the contacts is also obtained. 



Controversy on the observed electrical transport prop- 
erties of carbon nanotubes (CNT's) has been mostly due 
to our lack of control and understanding of their contact 
to the metallic electrodes. It has finally become clear that 
the contact influences critically the overall performance 
of the CNT and that it is crucial to lower the inherent 
contact resistance to achieve the definite understanding 
of the intrinsic electrical properties of CNT's [0, ^, In 
order to determine the relevant factors behind the con- 
tact resistance so that this can be pushed down to its al- 
leged quantum limit Rq = h/2e^ per CNT channel a big 
experimental effort has been made both in CNT growth 
and lithographic techniques [|, i |, 0, |, |, 0. While 
considerable progress in this direction has already been 
achieved, theoretical progress, on the other hand, lags 
behind in this important issue. 

The actual atomic structure of the electrode (and prob- 
ably that of the CNT) at the contact are unknown and, 
most likely, change from sample to sample when fabri- 
cated under the same conditions. Atomic-scale modeling, 
however, can still be of guidance to the interpretation of 
the experiments and to the future design of operational 
devices with CNT's. In this work we focus on the two 
key ingredients in this puzzle: The effect the atomic-scale 
geometry and the chemical nature of the electrode have 
on the transparency of the contact. We have studied 
open single- walled metallic (5,5) CNT's contacted in two 
representative forms (see Fig. |l|) to Al, Au, and Ti elec- 
trodes which are among the most commonly used metals 
in the experiments . From our ab-initio transport study 
we find that in CNT's contacted to Al and Au electrodes 
for end-contact geometry [see Fig. |l|(a)] the two CNT 
bands couple weakly to the electrodes. This allows us to 
resolve quasi-bound CNT states in the conductance and 
to estimate the magnitude of the degeneracy removal due 
to Coulomb blockade effects in a direct manner. More- 
over, we find that the two bands couple very differently 
to the electrodes (one of them is almost shut down for 




FIG. 1: The two contact geometries considered in this work: 
An open (5,5) carbon nanotube end-contacted to (111) sur- 
faces (a) and the same nanotube side-contacted (b). 



transport) and do not mix. For the side-contact geom- 
etry [see Fig. 0(b)] the coupling is the same for both 
bands, but similar in strength to the end-contact geome- 
try. Finally, our study presents the first direct numerical 
evidence of what has been hinted at on the basis of in- 
direct first-principles calculations [|ll[ O and what has 
recently been observed in experiments|7| : Early 3-d el- 
ements as Ti are probably the best choice for making 
high-transparency contacts to CNT's compared to more 
traditional metals such as Al and Au. Although perfect 
transparency at the contact is nerver achieved, our cal- 
culations indicate that properly engineered Ti contacts 
are a good bet for future perfect contacts to CNT's. 

From a theory point of view, the "contact" problem 
has been previously addressed [|l^, |l^, but only par- 
tially. The reason is that a full analysis of this problem 
requires the use of sophisticated state-of-the-art numeri- 
cal techniques to calculate electrical transport from first- 
principles 1 15, |l6|, |l^, where even the electrodes need to 
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be described down to the atomic level p^, |T^, ^ p2| . 
These techniques are currently under development. First 
of all, charge transfer at the contact, which aligns the 
chemical potentials of the electrodes and the CNT, needs 
to be evaluated self-consistently[^ Secondly, one 

needs to combine the ab-initio calculation with Lan- 
dauer's formalism!^. Recently, we have presented a 
very promising approach, termed Gaussian Embedded 
Cluster Method |2^, that allows us to address this 
problem in its full complexity. Our method is based on 
standard quantum chemistry calculations performed with 
the Gaussian98 code|2^]. A density functional (DF) cal- 
culation of a cluster comprising the CNT and a signif- 
icant part of the electrodes is performed (see Fig. [|). 
Next, the retarded(advanced) Green's functions associ- 
ated with the self-consistent hamiltonian or Fock opera- 
tor F of the cluster is modified to include the rest of the 
semi-infinite electrodes: 
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In this expression &^'> = + 1]^^'', where Sr(Sl) de- 
notes a self-energy operator that accounts for the part of 
the right (left) semi-infinite electrode that has not been 
included in the initial DF calculation [|2^, and / is the 
unity matrix. In a non-orthogonal basis, like those com- 
monly used in Gaussian98, the embedded cluster density 
matrix takes the form 
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where S is the overlap matrix, G^~^ is the retarded 
Green's function expressed in the non-orthogonal basis, 
and is the Fermi energy which is set by imposing over- 
all charge neutrality in the cluster. The density matrix is 
returned to Gaussian98 and the process is repeated until 
the procedure converges. The conductance can finally be 
calculated through the standard expression Pq]: 
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where Tr denotes the trace over all the orbitals in the 
cluster and where the matrices F^^ and are ^(S^^ — 

S^^) and ^(S^~^ — ^i^^), respectively. In order to single 
out the contribution of individual channels to the current 
one can diagonalize the transmission matrix T. 

Figure ||(a) shows Q around the Fermi energy for a (5,5) 
metallic CNT composed of = 10 carbon layers that has 
been end-contacted [Fig. 0(a)] to Al(lll) surfaces (the 
end-carbon-layer-surface distance has been optimized to 
a value of 1.8A)[|3^. Four resonances appear around the 
Fermi energy (set to zero). These resonances can be eas- 
ily traced back to four extended states of the isolated 
finite CNT|24|. Two of them {ki,k2) originate in the 
bonding (tt) band of the CNT and the other two (/cjf , /C2) 
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FIG. 2: (a) Conductance as a function of energy for an 
A = 10 (5,5) open metallic nanotube end-contacted to a 
Al(lll) surface [see Fig. ^(a)]. The nanotube-surface dis- 
tance has been optimized to a value of 1.8 Aand the Fermi 
energy is set to zero. Inset: Schematic band structure of the 
metallic nanotube showing the four states responsible for the 
resonances, (b) Transmission as a function of energy for the 
highest conducting channels. The symmetry of the two main 
channels is also shown. 



in the antibonding (tt*) band (see inset in Fig. ||). The 
resonances have different widths for different bands indi- 
cating that they couple very differently to the electrodes. 
Moreover, the two bands do not mix with each other. 
This is more clearly seen in Fig. ||(b) where we show the 
highest transmission eigenvalues of the transmission ma- 
trix. Two independent channels exhibit resonances in the 
energy window (~ 3.56^) around Ep where only the tt 
and TT* bands can contribute to transport. This result is 
consistent with the fact that tt* states, of large angular 
momentum, do not couple to the low-angular momentum 
states of the electrode, while tt states, of low angular mo- 
mentum, couple more easily [|l3|, Q. Notice that there 
is a charge transfer from the metal to the CNT, but this 
mainly localizes at the end carbon layer {^0.2 per carbon 
atom) and it does not affect the overall band positioning 
in the center of the CNT. 

The specific band assignment of the resonances is 
nicely confirmed by their evolution on the length of the 
CNT presented in Fig. ||. We have calculated the conduc- 
tance for AT = 8, 9, 10, 11, 12, and 13 carbon-layer CNT's. 
The opposite signs of the group velocity for the tt and tt* 
bands make the quasi-bound states belonging to the tt* 
band shift down in energies while those belonging to the tt 
band shift up as A" increases. As expected from a simple 
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FIG. 4: (a) Conductance as a function of energy for an = 
15 (5,5) open metallic nanotube side- contacted to a Al(lll) 
surface [see Fig. |l|(b)]. The nanotube-surface distance has 
been optimized to a value of 2.2 A. (b) Transmission as a 
function of energy for the three highest conducting channels. 
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FIG. 3: Conductance as a function of energy for an = 
8(a), N = 9(b), N = 10(c), N = 11(d), N = 12(e), and 
N = 13(f) (5,5) open metallic nanotube end-contacted to a 
Al(lll) surface [see Fig. 0(a)]. The Fermi energy has been 
set to zero. 



particle-in-a-box argument applied to finite CNT's[|4|, 
for N = 31^ where / is an integer, we should expect two 
states with the same wave vector kn but in different bands 
to coincide at the Fermi energy. Naively one should thus 
expect Q = 4e^//i[|^. Our results for the contacted 
= 9 and TV = 12 CNT's show otherwise: Two res- 
onances never coincide at the Fermi level. The reason is 
that Coloumb blockade prevents two (band and/or spin) 
degenerate quasibound states to be filled up at the same 
time and degeneracies are removed From Figs. ||(b) 
and (e) we estimate the charging energy to be ^ 0.3 eV 
in these CNT's which is smaller than the single-particle 
level spacing as confirmed by experiments 

If the interpretation of the different coupling strengths 
of the CNT bound states with the Al electrodes is cor- 
rect and angular momentum considerations are relevant, 
similar couplings should be expected for both bands if no 
axial symmetry is present. This is the case for the other 
contact geometry considered in this work [see Fig.|^(b)]. 
Figure ^ shows results for an = 15 CNT side-contacted 
to Al(lll) surfaces (the CNT-surface distance has been 
optimized to 2. 3 A). Conductance resonances come in 
pairs in the relevant energy window which is what is 



expected for an A^ = 15 CNT. More importantly, all 
of them present similar widths, confirming our expec- 
tations. Contrary to the previous geometry, localized 
end states |24| influence the coupling around leV for this 
contact geometry where mixing with the CNT extended 
states takes place. Our results for the coupling strength 
with Al contacts are consistent with previous studies 
where jellium models were considered as contacts [p!4|, 
and with those in Ref. 18, but we do not subscribe previ- 
ous ab-initio results presented in Ref. ^ based on what 
it seems to be more realistic contact models similar to 
ours. 

We now complete our study for end-contacted N = 10 
CNT's considering Au and Ti electrodes (see Figs. || and 
^). Several resonances are clearly visible close to the 
Fermi energy for the case of Au, but, in contrast to 
Al electrodes, it is difficult to identify specific extended 
states as we did above. This is in part due to the mix- 
ing of the TT and tt* bands with the end states which, 
in addition, induce extra channels in the conductance, 
although these channels are only relevant for transport 
in very short CNT's |3^. Apart from this, the coupling 
strength of the two bands is similar to that found for 
Al electrodes despite of the fact that the Mulliken popu- 
lation analysis reflects a minor charge transfer from the 
electrode to the CNT. In Fig. |^(b) we appreciate that 
the TT band coupling is also stronger than that of the tt* 
band. In contrast to Al and Au electrodes, where Q ex- 
hibits resonances, Q presents an oscillatory behavior for 
Ti around E^. This is accompanied, as the anticross- 
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FIG. 5: (a) Conductance as a function of energy for an = 
10 (5,5) open metallic nanotube end-contacted to a Au(lll) 
surface [see Fig. |l|(a)]. The nanotube-surface distance has 
been optimized to a value of 2.2 A. (b) Transmission as a 
function of energy for the highest conducting channels. 
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FIG. 6: (a) Conductance as a function of energy for an N = 
10 (5,5) open metallic nanotube end-contacted to a Ti(lll) 
surface [see Fig. |l|(a)]. The nanotube-surface distance has 
been optimized to a value of 1.8 A. (b) Transmission as a 
function of energy for the highest conducting channels. 



ings in the transmission eigenvalues reveal in Fig. ^(b), 
by band mixing. This result reflects, as suggested in Ref. 
|l|, that Ti couples differently to the CNT (due to the 
presence of (i-states at the Fermi energy) and forms a 
better contact (the charge transfer is ^ 0.4 electrons per 



C atom at the end layer). At this point, however, we can 
only speculate on the possibility of perfect transparency 
for other Ti electrode geometries. 
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